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^ ■ ABSTRACT 
O 

We investigate electromagnetic buoyancy instabihties of the electron-ion plasma with 
the heat flux based on not the magnetohydrodynamic (MHD) equations, but using the 
multicomponent plasma approach when the momentum equations are solved for each 
species. We consider a geometry in which the background magnetic field, gravity, and 
stratification are directed along one axis. The nonzero background electron thermal flux is 
taken into account. Collisions between electrons and ions are included in the momentum 
equations. No simplifications usual for the one-fluid MHD-approach in studying these 
instabilities are used. We derive a simple dispersion relation, which shows that the thermal 
flux perturbation generally stabilizes an instability for the geometry under consideration. 
This result contradicts to conclusion obtained in the MHD-approach. We show that the 



- 2 - 



reason of this contradiction is the simphfied assumptions used in the MHD analysis of 
buoyancy instabihties and the role of the longitudinal electric field perturbation which is 
not captured by the ideal MHD equations. Our dispersion relation also shows that the 
medium with the electron thermal flux can be unstable, if the temperature gradients of ions 
and electrons have the opposite signs. The results obtained can be apphed to the weakly 
coUisional magnetized plasma objects in laboratory and astrophysics. 

Keywords convection - instabilities - magnetic fields - plasmas - waves 
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1. INTRODUCTION 

Thermal effects resulting in instabilities, transport, heating, structures forming, and so 
on play an important role in dynamics of different plasma objects in laboratory, space, and 
astrophysics. For example, the ion-temperature-gradient-driven modes (Kadomtsev and 
Pogutse 1965) are used for explaining anomalous transport in tokamak plasma experiments 
(Dimits et al. 2000; Garbet 2001). Thermal conductivity influences on the Rayleigh-Taylor 
instability in inertial fusion (Betti et al. 1998; Canaud et al. 2004; Lindl et al. 2004), 
on the surface of the Sun (Isobe et al. 2005), in supernova (Fryxell et al. 1991), and 
other astrophysical objects. Thermally stratified fluids can be buoyantly unstable in the 
gravitational field. In astrophysics, this process may, for example, operates in the stellar 
interiors (Schwarzschild 1958), accretion disks (Balbus 2000, 2001), neutron stars (Chang 
and Quataert 2010), hot accretion flows (Narayan et al. 2000, 2002), galaxy clusters, 
and intracluster medium (ICM) (Quataert 2008; Parrish et al. 2009; Ren et al. 2009, 
2010a). Analogous instabilities also exist in the neutral atmosphere of the Earth and ocean 
(Gossard and Hooke 1975; Pedlosky 1982). Diversity of environments in which buoyancy (or 
convective) instabilities may have the significant role, leading to turbulence and anomalous 
energy and matter transport, makes these instabilities an important object for analytical 
and numerical explorations. 

The crucial role of convection in the transport of energy, for example, in stellar 
interiors is a well-known physical process (Schwarzschild 1958). However, theoretical efforts 
to understand convective energy transport in the dilute and hot plasmas such as galaxies 
clusters and ICM (Sarazin 1988) have lead to some results over recent years. As it is 
known, majority of the mass of a cluster of galaxies is in the dark matter. However, around 
1/6 of its mass consists of a hot, magnetized, and low density plasma known as ICM. 
The electron number density is rig ~ 10~^ to 10"^ cm~^ at the central parts of ICM. The 
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electron temperature Tg is measured of the order of 1 — 15 keV, though the ion temperature 
Tj has not yet been measured directly (Fabian et al. 2006; Sanders et al. 2010). The 
magnetic field in ICM is estimated to be in the range 0.1 — 10 /iG depending on where the 
measurement is made (Carilli and Taylor 2002). This implies a dynamically weak magnetic 
field with /3 — SnUeTe/ B'^ 200 — 2000. Under conditions given above, the Larmor radius 
of electrons and ions (Tj ~ Tg) is many orders of magnitude smaller than the mean free 
path. Thus, the ICM is classified as a weakly coUisional plasma (Carilli and Taylor 2002) 
possessing anisotropic transport due to the magnetic field. 

In recent past, for such plasmas in the framework of the ideal MHD supplemented by 
an anisotropic heat flux along the magnetic field, there were found some new convective 
instabilities for the case when a heat flux plays the significant role (Balbus 2000, 2001; 
Quataert 2008; Ren et al. 2009, 2010a, 2010b). One of these instabilities, at the absence of 
the background thermal flux, has been shown to arise when the temperature increases in 
the direction of gravity which is perpendicular to the background magnetic field. This is 
so-called the magnetothermal instability (MTl) (Balbus 2000, 2001). The other instability, 
the heat buoyancy instability (HBI) (Quataert 2008), has been found to arise at the 
presence of the background heat flux when the temperature decreases along gravity parallel 
to the magnetic fleld. Anisotropic dissipative effects have been included by Ren et al. 
(2010a, 2010b). 

Theoretical models applied for study of buoyancy instabilities in astrophysical objects 
with a heat flux are based on the one-fluid ideal (Balbus 2000, 2001; Quataert 2008; Ren 
et al. 2009; Chang and Quataert 2010) and nonideal (Ren et al. 2010a, 2010b) MHD 
equations. By using of these equations one can comparatively easily to consider any 
problems. However, the ideal MHD docs not capture some important cfi^ccts which can 
be taken into account by using a multi-fluid plasma approach. One of such effects is the 
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nonzero longitudinal electric field perturbation along the background magnetic field. An 
importance of involving this component due to multi-fluid effects and shortcomings of the 
ideal MHD were emphasized, for example, for the acceleration of solar flare electrons by 
inertial Alfven waves (McClements and Fletcher 2009), at consideration of structures of 
electromagnetic fields and plasma flows in pulsar magnetosphere (Kojima and Oogi 2009), 
for the acceleration of relativistic ions, electrons, and positrons in shock waves (Takahashi 
et al. 2009), and in a gyrofluid description of Alfvenic turbulence (Bian and Kontar 2010). 
As we show here, the contribution of currents due to this (small in the present case) 
parallel electric field to the dispersion relation can be of the same order of magnitude as 
that due to transverse electric field components. Besides, the MHD equations do not take 
into account the existence of various charged and neutral species with different masses and 
electric charges and their collisions between each others and therefore can not be applied 
to multicomponent systems. In some cases, the standard methods used in the MHD lead 
to conclusions that are different from those obtained by the method using the electric field 
perturbations (the E-approach). One such an example concerning the contribution of the 
electron-ion collisions to the dispersion relation for the MHD waves in the two-component 
magnetized plasma was considered by Nekrasov (2009c). A multicomponent approach has 
been used in (Nekrasov 2008, 2009a, 2009b, 2009c, 2009d), where the streaming instabilities 
of rotating astrophysical objects (accretion disks, molecular clouds and so on) have been 
investigated. 

A study of buoyancy instabilities with the electron heat flux by the multicomponent 
E-approach has been performed by Nekrasov and Shadmehri (2010). The geometry has 
been considered in which the gravity is perpendicular to the background magnetic field and 
the background heat flux is absent. Solution of the dispersion relation obtained in this 
paper differs from solution of the same problem found from the ideal MHD (Balbus 2000). 
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In this paper, we apply a multicomponent approach to study buoyancy instabihties in 
magnetized electron-ion plasmas with the background electron thermal flux. We consider 
the geometry in which the gravity, stratification, background magnetic field and thermal 
flux are all directed along one (z-) axis. For generality, we include collisions between 
electrons and ions in the momentum equations. At the consideration of the perturbed heat 
flux, we adopt that cyclotron frequencies of species arc much larger than their collision 
frequencies. Such conditions are typical for many laboratory, space, and astrophysical 
plasmas. In this case, the heat flux is anisotropic and directed along the magnetic fleld 
fines (Braginskii 1965). However in other respects, the relation between the cyclotron and 
collision frequencies is arbitrary in the general expressions for the perturbed values. We 
derive the dispersion relation for cases, in which the background heat flux is present or 
absent. This gives a possibility to compare these two cases. Solutions of the dispersion 
relation are discussed. 

The paper is organized as follows. In Sect. 2, the fundamental equations are given. 
An equilibrium state is considered in Sect. 3. Perturbed ion velocity, number density, and 
thermal pressure are obtained in Sect. 4. In Sect. 5, we consider the perturbed velocity 
and temperature for electrons. Components of the dielectric permeability tensor are found 
in Sect. 6. Dispersion relation is derived and considered in the coUisionless and colfisional 
cases in Sect. 7. Discussion of the results obtained and comparison with the MHD results 
are provided in Sect. 8. In Sect. 9, we give conclusive remarks. 



2. BASIC EQUATIONS 



We start with the following equations for ions: 




(1) 
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the momentum equation, 

_^ + V • n,v, = 0, (2) 

the continuity equation, and 

^ + • Vpi + 7PiV • = 0, (3) 
the pressure equation. The corresponding equations for electrons are: 

= -Z^+g^E+^Ve X B-mei/ei(Ve-Vi), (4) 

ne c 

^ + V-neVe = 0, (5) 

^ + Ve • Vpe + 7PeV • Ve = A - (7 - 1) V • q„ (6) 

c)T A X 

+ Ve • VTe + (7 - 1) TeV • Ve = - - (7 - 1) "V • q„ (7) 

at He rie 

the temperature equation, where Qg is the electron heat flux (Braginskii 1965). We neglect 
inertia of the electrons. In (l)-(7), qj and rrij arc the charge and mass of species j = i,e, 
Vj is the hydrodynamic velocity, rij is the number density, pj — UjTj is the thermal 
pressure, Tj is the temperature, i/jg (j^ei) is the collision frequency of ions (electrons) with 
electrons (ions), g is gravity, E and B are the electric and magnetic fields, c is the speed 
of light in vacuum, and 7 is the adiabatic constant. At the consideration of the electron 
heat flux, we will assume the electrons to be magnetized when their cyclotron frequency 
^ce — QeB/rrieC ^ J^eei^ei), whcrc ^eei^ei) IS the clectron-clectron (ion) collision frequency. 
In this case, the electron thermal flux is mainly directed along the magnetic fleld, 

qe = -Xeb(b-V)Te, (8) 

where Xe is the electron thermal conductivity coefficient and h = 'B/B is the unit vector 
along the magnetic field (Braginskii 1965). However in the momentum equations (1) and 
(4), we keep for generality an arbitrary relation between ujci{ujce) and I'lei'i^ei) i^d — QiB/rriiC 
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and fie are the ion cyclotron and ion-electron collision frequencies, respectively), having 
in mind that some expressions obtained below can be applied for coUisional objects. The 
term A compensates the temperature change as a result of the equilibrium heat flux. We 
take only into account the electron thermal conductivity by (8), because the corresponding 
ion conductivity is considerably smaller (Braginskii 1965). For generality, we assume the 
electron and ion temperatures to be different. However, we do not involve, for simplicity, 
the terms describing the energy exchange between ions and electrons in (3), (6), and (7). 
Thus, our treatment is available for cases in which such an exchange is not effective or when 
there is a strong temperature coupling of species. In the last case, one can set Tj ~ T^. This 
issue is considered in more detail in Sect. 8. 

Electromagnetic equations are Faraday's law 

V X E = 

and Ampere's law 

V X B 

where j = qjUjVj. We consider the wave processes with typical timescales much larger 
than the time the light spends to cover the wavelength of perturbations. In this case, 
one can neglect the displacement current in (10) that results in quasineutrality both in 
electromagnetic and purely electrostatic perturbations. The magnetic field B includes the 
background magnetic field Bq, the magnetic field ^ocur of the background current (when it 
presents), and the perturbed magnetic field. 

3. EQUILIBRIUM STATE 

At first, we consider an equilibrium state. We assume that background velocities are 
absent. In this paper, we study configuration in which the background magnetic field. 



IdB 
c~dt 



=— J, 



(9) 



(10) 
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gravity, and stratification are directed along the 2;-axis. Let, for definiteness, g be g = — 
where g' > and z is the unit vector along the z-direction. Then, (1) and (4) give 

1 dpio Qi 

Qi = 7^ = 9 Eo, (11) 

miTiio oz rrii 

Qe = ^ = — ^0, (12 

miTieo OZ rrii 

where (and below) the index denotes equilibrium values. Here and below, we assume that 
Qi = —Qe- For convenience of notations, we do not use that n^o = for the two-component 
plasma up to a point where it will be necessary. We see that equilibrium distributions of 
ions and electrons influence on each other through the background electric field Eq. In the 
case riio = rieo and T^o = Tgo, we have Qi = Qe = g/'^- Thus, wc obtain Eq = mig/2qi. The 
presence of the third component, for example, of the cold dust grains with the charge qa and 
mass nid ^ rrii results in other value of £'0 = Tnag/qd- In this case, the ions and electrons 
are in equilibrium under the action of the thermal pressure and equilibrium electric field, 
being Qi ~ -g^, if qiUid > qdrrii. 



4. LINEAR ION PERTURBATIONS 



Let us write (l)-(3) for ions in the linear approximation, 

— 1 hijiH Vji X Bo, (id) 

dt miTiio miriionio rriiC 

dna duio , . 

-^ + Viu-^ + nio\/ -Yii^O, (14) 

dpii dpio , . 

-Q^ + Viu-Q^ + 'jPioV ■Vii = 0, (15) 



where 



Fa = ^Ei - (va - Vei) , (16) 



and the index 1 denotes the perturbed variables. Below, we solve these equations to find 
the perturbed velocity of ions in an inhomogeneous medium. 
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4.1. Perturbed velocity of ions 

Applying the operator d/dt to (13) and using (14) and (15), we obtain equation 
containing only the ion velocity 

^= -9iVv,,, + ^ [(7 - 1) (Vp,o) + 7P.oV] V • va + ^ + X Bo- (17) 

at^ miTiiQ at rriiC at 

Prom this equation, we can find solutions for the components of Vji- For simplicity, we 

assume that d/dx = because a system is symmetric in the transverse direction relative to 
the 2;- axis. The x-component of (17) has the simple form 

—^—Fiix + IjJciVily (18) 

Here uici — qiBo/rriiC. For the y-component of (17), we obtain: 

where, Csi = {'yTio/rrii)^^'^ is the ion sound velocity. Using (18), equation for vny is given 
from (19) as follows 

2\ ^ dPa 

Q^+^cij Vily - Qily^-g^. (20) 



Then from (18), we obtain 



d 



Uridt 



92 2 , 



dy 



(21) 



In (20) and (21), we have introduced the following notations: 

Pii = -QiViiz + 4V • va, (22) 

Qilx = ^ciFily H '' ^^"^-^ 

Qily = —l^ciFilx H (24) 
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The value Pa defines the ion pressure perturbation (see 15). We see from (20) and (21) 
that the thermal pressure effect on the velocity Vn^ is much larger than that on Vny when 
d/dt <^ ujci. 



The z-component of (17) takes the form 



d f dviu 



dt\dt 



oz 



„2 d 



V • va- 



(25) 



To obtain equation only for Viu, we need to express V • Va through vnz- Differentiating 
(20) with respect to y and using expression (22), we find 

dQily 



where the following operators are introduced: 



^2 



dy 



(26) 



(27) 



(28) 



We can derive equation for the longitudinal velocity Vjiz, substituting V • Vji found from 
(26) into (25), 



, , , dQ 

J^3Viiz—J^l—Q^ r -L4— ^— , 



(29) 



where operators L3 and L4 have the form 



Qy2 dz"^ J dt"^ 



2 2 ^ 



(30) 



q2 \ d dL 

+ 79. 1 + i^l, j 9j + + (1 - 7) al-g-i, 



4(|_||). (31) 

For obtaining expression (30), we have used expressions (27) and (28). 
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It is easy to see that at the absence of the background magnetic field and without 
taking into account electromagnetic perturbations and collisions with electrons (the right 
hand-side of 29), equation L^Vn^ = describes the ion sound and internal gravity waves. In 
this case, a sum of the last two terms on the right hand-side of expression (30) is equal to 
— CgjCu^j^, where Ubi is the (ion) Brunt- Vaisala frequency equal to 



Thus, we have obtained a result corresponding to perturbations in the neutral atmosphere 
(Gossard and Hooke 1975). However, we see that the existence of the background magnetic 
field considerably modifies the operator L3. We note that the right hand-side of (29) 
describes a connection between ions and electrons through the electric field Ei and collisions. 



So far, we did not make any simplifications and all the equations and expressions 
obtained above are given in their general form. This permits us to investigate different 
limiting cases. Further, we consider perturbations with a frequency much lower than the 
ion cyclotron frequency and the transverse wavelengths much larger than the ion Larmor 
radius. Such perturbations are of interest for both laboratory and astrophysical plasmas. 
Besides, we investigate a part of the frequency spectrum in the region lower than the ion 
sound frequency. Thus, we set 



^li = ^ (7 - 1) 5i + 




(32) 



dz 



4.2. Specific case for ions 




• ^ 



(33) 
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In this case, operators (27), (28), (30), and (31) take the form 



T 2 T 2 9 



(34) 



d 



L^ = (l-7)c/i + 4 



d 



d_ 
dz 



dz 



Also, an additional condition 



2 2 dc% 

LU ■ ^ C • — 

"dt^ dz dy-'dz 



(35) 



must be satisfied for operator L3 to have a given form (34). The small corrections in 
operators L3 and L4 are needed to be kept because some main terms in expressions for the 
ion and electron velocities are equal to each other (see below). Therefore, when calculating 
the electric current, these main terms will be canceled and small corrections to velocities 
will only contribute to the current. The applicability of condition (35) and other conditions 
used below will be discussed in Sect. 8. 

For cases represented by inequalities (33) and (35) when the operators Lj, i = 1,2,3,4, 
have the form (34), equations for vn^ and V • Vji become 



d 



d 



""''dz ^^V dz dt'' 



dR 



ilz 



dt 



d_ 
dz 



dQ 



ily 



dz ujiidy 



(36) 
(37) 



4.3. Ion perturbations in the Fourier transformation 

Calculations show that some main terms in expressions for (when calculating 
the current), V ■ Vji and Pn are canceled. Therefore, the small terms proportional to 
inhomogeneity must be taken into account. To do this correctly, we can not apply the 
Fourier transformation to (36) and (37) to find the variable Pn. However, firstly, we should 
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apply the operator d/dz to this variable for using (36). It is analogous to obtaining the 
term dc^g/dz in expression (32) for the Brunt- Vaisala frequency. After that, we can apply 
the Fourier transformation in a local approximation assuming the linear perturbations to 
be proportional to exp(ikr— ia;i). As a result, we obtain for the Fourier-components vnzki 
k • Vjifc, and Pjife, where k — (k,a;), the following expressions: 

k • vafe= -ly-j- 1 - ^7— 2- Fiizk + ^7r;j-^Q«i2/fc' (39) 

^zCgi \ ^z^si/ ^zCsi^ci 



, - iLz^ - , 2„2 



(7 - 1) + 



dc 



dz 



(7-l)^. + ^-cu^^ 



i^azife (40) 

Qilyk-i 



where gi/kzC^^ <^ 1. In expressions (38) and (39), we have omitted additional small terms 
Qiiyk which are needed for calculation of Puk. 

When calculating the current along the z-axis, the main term ~ Qnyk in (38) will be 
canceled with the corresponding electron term. The contribution of the first term ~ Fn^k 
to this current has, as we will see below, the same order of magnitude for the buoyancy 
instabilities as contribution of the term ~ giQnyk, i-e. -Piizfc ~ {^y9i/^^ci) Qnyk- The same 
relates to expressions (39) and (40). Thus, the longitudinal electric field perturbation Eiz 
containing in Fn^ must be taken into account. However, in the ideal MHD, this field is 
absent. We see from expressions (38) and (39) that V • Vji ~ {gi/ c^g^Vnz- This relation is 
the same as that for internal gravity waves in the Earth's atmosphere (see, e.g., Nekrasov 
1994). Using expression (40), we can find velocities v^yk and Vnxk from (20) and (21), 
corresp ondingly. 
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4.4. Perturbed ion number density and pressure 



It is followed from above that V • Va ~ Vuz/H, where H is the inhomogeneity 
scale height {H ~ c^si/dd- Thus, the last two terms in (14) and (15) are of the same 
order of magnitude. Let us find the perturbed ion number density and pressure in the 
Fourier- representation. Prom (14), (38) and (39), we obtain 



nuk 
riio 



Fa 



zk 



Q 



ilyk- 



Equation (15) gives dpn/dt — —minioPn. Thus, we have, using (40), 



Piik 

PiO 



7 



Filzk + 



(7-1).. + ^ 



to 



lyk- 



(41) 



(42) 



Comparing (41) and (42), we see that the relative perturbation of the pressure due 
to the transverse electric force Qnyk is much smaller than the relative perturbation of 
the number density. However, these relative perturbations as a result of the action of the 
longitudinal electric force F^zk have the same order of magnitude (see Sect. 4.3). Thus, 
Piik/Pio ~ "nnk/niCj. This result contradicts a supposition Puk/Pio <^ nuk/riio adopted in the 
MHD analysis of buoyancy instabilities (Balbus 2000, 2001; Quataert 2008). From results 
obtained below, it is followed that, as we already have noted above, the both terms on the 
right hand-side of (41) have the same order of magnitude. 



LINEAR ELECTRON PERTURBATIONS 



Equations for the electrons in the linear approximation are the following: 



= \ h Fel H Vel X Bo, (43) 

UeO rieo rieo C 
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-^ + VeU-gf + neoV-Vel^O, (44) 

+ '"eiz-Q^ + TPeoV • Vei = - (7 - 1) V • q^i, (45) 

^ + ^eU^ + (7 - 1) TeoV • Vel = - (7 - 1) —V • q.^, (46) 

qei = ~"i^eo-^ - boXeo-^ - hoXel-Q^, (47) 

Fei = QeEi - nieUei (Vel " Va) ■ (48) 

In (47), Xei = 5XeoTei/2TeQ is the perturbation of the thermal flux conductivity coefficient 
Xe which is proportional to Te^^ (Spitzer 1962; Braginskii 1965). The perturbation of the 
unit magnetic vector bi is equal to 61^ ^ = Bi^^y/Bo and bu = 0. The thermal flux in 
equihbrium is qeo = -boXeo^- 

We have seen above at consideration of the ion perturbations that the terms ~ 1/if^ 
are needed to be kept (see the last term in 40). Therefore, we also keep such terms for the 
electrons. 



5.1. Equation for the electron temperature perturbation 

Let us now find equation for the electron temperature perturbation. Expression V • q^i, 
where qei is defined by (47), is given by 

'^^ dy dz ^° dz Bq dy ^° dz^ dz dz dz'^ 

Substituting this expression into (46), we obtain 

rp 9TeQ . , XeodTeodBiy 

D,T^, = ---^ - (7 - 1) TeoV • ve, + (7 - 1) (50) 



where the operator Di is defined by 



dt Ueo V dz^ dz dz dz'^ 



(51) 
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5.2. Perturbed velocity and temperature of electrons 

We further find equations for components of the perturbed velocity of electrons. The 
x-component of (43) has a simple form, i.e. 

Vely = Felx, (52) 

rrieUce 

where oUce — qeBo/rrieC. Applying the operator d/dt to the y-component of (43) and using 
(45) and (49), we obtain 

dt \ ^^"^ rrieUJce ""^^ ) ^ci dy me^ce^eo dz Body^ 



rUeUJce V dt J dy 
Here 

Pel = -QeVelz + C^V • Vgl, (54) 

where c^^ = jpeo/ ^ineo- The variable Pgi is analogous to Pn (see 22) and defines the 
electron pressure perturbation. But for electrons, their pressure perturbation is also affected 
by the thermal conductivity (see 45). The 2;-component of (43) takes the form 

1 dpel 1 dpeOUel , . 

0= ^ \ ^ h Feu. (55) 

UeO OZ UeO OZ 71^0 

We can express V • Vgi which is contained in (54) through Veu, using (52), 

^ dz meUJce dy 
This expression for electrons is analogous to that for ions (see 37) . 

We further consider perturbations with the dynamical frequency d/dt satisfying the 
following conditions: 

> — > J_^X^_^ 

Ueo dz'^ dt dz dz 

The first inequality (57) means that the thermal conductivity is the dominant mode of the 
thermal transport (Balbus 2000; Quataert 2008). Under the second condition (57), we can 
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neglect the inhomogeneity and perturbation of the thermal flux conductivity coefficient 
in the temperature equation (50) (see 47 and 51). Obviously, the term proportional to 
d'^Xeo/dz^ in expression (51) can also be neglected. In this case, the small correction 
proportional to d/ dt in the temperature equation (50) which will be necessary for calculation 
of the electron velocity (see below) will be larger than that ~ dxeo/dz. We further apply 
the operator d/dt to (55) and use (44), (45), (49), and (56). As a result, we obtain 



d 



dz 



rriidt 
+ (7-1) 



+ 



(l-7)^?e + C^ 



d_ 



XeO 



dT,o 1 d^B.y 



1 dF,,, 

nieUce dy 



(5^ 



miUeo V -^0 dydz ' dz^ 
Equation for the temperature perturbation under conditions (57) has the form 



(7-1) 



Xeod^ _ d_ 
UeQ dz^ dt 



Tel = v^u^ + (7 - 1) Teo - (59) 

oz \ oz rUeUJce oy ' 

I 1 X XeO ^T'eo dBxy 
- 7-1 r: , 

Ueo OZ Body 



where we have used (56). 



To find equation for Vgi^, we substitute T^i from (59) into (58). Taking into account a 
contribution of the term dTgi/dt and carrying out some transformations, we obtain 



f dY'd'^F,^, , 1 d^F,^, 

+ a„.a^2 (60) 



dz^ Teodzdt XeO \dz J T^odf^ rUeOJce dydz"^ 

hy 



1 ( , dcl\ 1 a^Fei, aTeO 1 a^Si 



^ \ ^ dz J nieUJce dydz Teodz Bo dydt 

The correction proportional to dFeix/dt is absent. The last term on the right hand-side of 
(60) is connected with the background electron thermal flux. 
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Prom (59) and (60), we can find equation for the temperature perturbation 



Ueo dz \ dz'^ dz Body 



eO 



/ . N dcl„ 

(7-l).e+^ 



(61) 



9^ 

7^ I TT 



eO 



12/ 



9^ V^-^/ Bodydt 

It is followed from results obtained below that all terms on the right-hand side of (61) 
(except the correction ~ d'^F^i^/ dt^) have the same order of magnitude (see Sect. 4.3). 
The left-hand side of this equation is larger (see conditions 57). Thus, the temperature 
perturbation in the zero order of magnitude can be found by equaling the left part of (61) 
to zero. However, the right part of this equation is necessary for finding the transverse 
velocity perturbation v^ix (see below). 

To find the velocity Weia;, we need to calculate the value P^i (see 53 and 54). Performing 
calculations in the same way as that for ions (see Sect. 4.3), we obtain 



dz'^ 



d_ 

'Yz 



+ (7 - 1) 5e + 



dz 



dFei, , dVei 



dz 



rriidt 



+ 



(62) 



+ 9e 



dc^ ^ 



dF, 



elx 



rUeOJce dy 



where we have introduced notation connected with the thermal conductivity. 



Ki = (7-1) 



XeO 



dT,o 1 dB^y 



+ 



rriineo \ dz Bq dy ' dz"^ 
Equation (62) can be re-written in the form which is convenient for finding the velocity 
Veix- Using (61), we obtain 



(63) 
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(64) 
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elx 
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eO 



Teodz V dz 



d 



d^B 



Bodydt 



It is easy to see that (53) has the form 




Fely) — TT- {Pel — Kl) ■ 



(65) 



Thus, the main contribution of the flux described by (63) does not influence on the electron 
dynamics. Applying operator d^/dz^ to (65) and using (64), we find an equation for the 
perturbed velocity Veix- 

6. FOURIER CURRENT COMPONENTS 

6.1. Fourier velocity components of ions and electrons 

Let us give velocities of ions and electrons in the Fourier- representation. From (20), 
(21), and (40), we have 



'^ilyk— 2 ^ U2, ,2 Qilyk 

^ci I l^z^ci 

Here and below, we have introduced notations 




(67) 



(66) 



c: 



.2 

'si,e . 



1 



(7 - 1) 9i,e + 




(68) 



The velocity vnzk is given by (38). 



For electrons, using (64) and (65), we find 




(69) 



where the following notation is introduced: 




(70) 
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Equation (60) gives us 

Velzk — -J ^elxk — ^To elxk + — — — {( I) 

- i ^ (l + icu^^^ F 
The velocity v^iy is defined by (52). 



6.2. Fourier electron velocity components at the absence of the heat flux 

To elucidate the role of the electron thermal flux, we will also consider the dispersion 
relation when this flux is absent. Therefore, we also give here the corresponding electron 
velocity components: 

_ ^y^eQe 11 ky ( aA 1 

k^LOUcimeUJce rUleUJce K\ k^ J nieUJce 

Veizk^-r 1 - F^ixk-i'-pr^ 1 - — 5- -feUfe- (73) 

Comparing expressions (69) and (71) with these equations, we see that the thermal flux 
essentially modifies the small terms in the electron velocity under conditions (57). 



6.3. Fourier components of current 

We find now the Fourier components of the linear current ji = qiUioVn + ge^eoVei- It is 
convenient to consider the value Am^i/u. Using expressions (38), (52), (66), (67), (69), and 
(71), we obtain the following current components: 

47ri 



-jlxk — O^xxEixk + iO-xyEiyk — O^xzElzk (74) 
^xx (^ilxfe ^elxfe) ibxy {Vilyk ^elyfe) ~l~ ^xz (^ilzfc ^elzfe) ) 
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47ri 

jlyk — ~iO'yxEixk + O'yyEiyk — dyzEizk (75) 

+ ibyx {Vilxk — Vgixk) — byy (wjiyfc " Wgli/fc) + byz {vnzk ^ ^eUk) ) 

47ri 



a; 



-jlzk — —O'zxElxk ~ (^zyElyk + dzzEiz (76) 
+ fe^x {Vilx - Velx) + bzy {Viiy - Vely) - bzz {vhz - Velz) ■ 



When obtaining (74)-(76), we have used notations (16), (23), (24), and (48) and equahties 
Qe — —Qi, rieo — JT-iOj 'iTT'e^ei — 'fni^ie- We also have substituted Biy). by {kzc/uj)Eixk, using 
(9). The following notations are introduced in (74)- (76): 

"""" ~ ul kl V k^ k^ nodz J' ^ ^ 
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^ci ^ci OJOJci Kz \ Csi ^ eaOZ 



i-2 



z \ se SI 



and 



■'XX 



^pi^ie mi k'^ 


0- 


ul q, kl 


< ky A 
uucikl \ 


_ El 


^si. 


rrii 









kl Qitti + aecl^be 
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(78) 



bzx = —^'-1 1^,-^1 —u 



le-i 



"13 

Here ujpi = {Annioqf /mi)^^'^ is the ion plasma frequency and k'^ = k^ + kl- The terms 
proportional to T*q are connected with the background electron thermal flux. 

Calculations show that to obtain expressions for a^j without thermal flux, using 
electron velocities (72) and (73), we must change b^ by Qe/c^sei P^^ ^e*o ~ ^1 ^Skkis 7=1 
in terms axz and a^^. 
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6.4. Simplification of collision contribution 

An assumption that electrons are magnetized has only been involved by neglecting the 
transverse electron thermal flux. In other respects, a relationship between uice and v^i or 
ujci and Vie (that is the same) can be arbitrary in (74)- (76). We further proceed by taking 
into account that ui <^ujci- In this case, we can neglect the coUisional terms proportional to 
hxy and hyx- However, the system of equations (74)-(76) stays sufficiently complex to flnd 
ji through El. Therefore, we further consider the case in which the frequency cj and wave 
numbers satisfy the following conditions: 

Vl F » ^ KH- ' ^''^ 



where 



2 ^ ^si^se /gQN 

It is clear that conditions (79) can easily be realized. In this case, the current components 
are equal to 

4:7Tt 

jlxk ^xx-^lxk ~l~ ^^xy-^lyk ^xz-^lzki (^l) 

(jO 

Am . •IT' ZT- ZT- 

Jlyk — '^^yx-t^lxk ~l~ ^yy-^lyk ^yz^lzki 

(jj 

3\zk ^zx-^lxk ^zy-^lyk ~l~ ^zz-^lz- 

OJ 

Components of the dielectric permeability tensor £y are the following: 

_ . Vie ky {aj - 7ae) _ Vje ky {aj - ^ae) , . 

Sxx - a,, + i^^^ ^2 (1 _ ^^^) ^-y - + ^2 (1 _ ^^^) 

'T^yx — ^yx ~~^Tr7\ ■ J \ ' ^yy ~ ^yy> 



CI 

Qjyz ^zx ^'zy (i'-'7 



- (1 - z4) ' (1 - z4) ' (1 - Id,) ' (1 - z4) ^ 

where we have used notations (78). Parameter d,. 
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defines the coUisionless, <^ 1, and coUisional, dz ^ 1, regimes. Below, we derive the 
dispersion relation. 

7. DISPERSION RELATION 

Using (9) and (10) in the Fourier-representation and a system of equations (81), we 
obtain the following equations for the electric field components: 

(^^ — ^xx) Eixk — iSxyEiyk + £xzEizk = 0, (84) 
iSyxEixk + {n^z ~ ^yy) E\yk + {—nyUz + £yz) Eizk — 0, 
£zxEixk + {—nyUz + £zy) Eiyk + {riy — Ezz) E^k = 0, 

where n — 'kc/cu. The dispersion relation can be found by setting the determinant of the 
system (84) equal to zero. In our case, the terms proportional to Sxy and Syx can be 
neglected. As a result, we have 

(n £xx) Ij^y^yy ~^ (j^z ^yy} ^zz '^y'^z i^yz ~l~ ^zj/) ~^ ^yz^zy^ 

+ {nl - Syy) ExzSzx = 0. (85) 

This dispersion relation can be studied for different cases. In subsequent sections, we 
consider both the coUisionless and coUisional cases. 

7.1. CoUisionless case 

We assume that condition 



(86) 
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is satisfied. Then, using notations (77) and (82), we reduce the dispersion relation (85) to 
the form 



where ca — Bq/ [AT^minioY/'^ is the Alfven velocity and 



(87) 



Q = Qitti + C^^ttehe + C^^aeTz— - + (Cj - 7ae) 6e - — + T^T^ 



(88) 



When obtaining (87), we have used the condition kyCyoJ^^ <S 1. We see that there are 
two wave modes. The first mode, u"^ = k'^c\, is the Alfven wave with a polarization of 
the electric field mainly along the y-axis (remind that the wave vector k is situated in the 
y — z plane). This wave does not feel the inhomogeneity of medium. The second wave has 
a polarization of the magnetosonic wave, i.e. its electric field is directed mainly along the 
X-axis (see below). This wave is undergone by the action of the medium inhomogeneity 
effect. The corresponding dispersion relation is 



(89) 



Expression (88) can further be simplified, using (11), (12), (68), (70), and (80). As a result, 
we obtain 



7 



(74 + cL) m'f 



(7 - 1) niig + 7 



d{Tio + T,o) 
dz 



mi9 + 



dz 



(90) 



We have pointed out at the end of Sect. 6.3 what changes must be done in expressions 
(77) and (78) to consider the case without the electron heat flux. This case follows from 
(90), if we omit the term d (Teo + T*q) /dz and put 7 = 1 in the first multiplier. Then 
becomes (Q^ Qf) 



(7-1)5 + 



^(4 + 4) 

dz 



(91) 
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This is the Brunt- Vaisala frequency. Comparing expressions (90) and (91), we see that 
the heat flux stabihzes the unstable {fl'l < 0) stratiflcation. We also see from (90) that 
the background heat flux (~ T*q) has no a fundamental importance. If the temperature 
decreases in the direction of gravity {dTi^eo/dz > 0), the medium is stable. Solution (90) 
describes an instability only when 

7-1 ^ dTo ^1 
— — miQ < — — < -miQ, 
27 ^ dz 2 ^' 

where Tjo ~ Teo = Tq. We also note that can be negative if gradients of Tjo and Teo have 
different signs. 

The dispersion relation (87) with deflned by (90) considerably differs from the 
dispersion relation obtained in the framework of the ideal MHD (Quataert 2008). The 
reasons of this are discussed in Sect. 8. 



7.2. Collisional case 

We proceed with the collisional case when 

^^^^ 

In this limiting case, the dispersion relation takes the form — s^x ~ and we obtain 
again (89). Thus, the dispersion relation is the same for both the coUisionless and coUisional 
cases. We note that this result has also been obtained for the case in which gravity is 
perpendicular to the magnetic fleld (Nekrasov and Shadmehri 2010). 
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7.3. Polarization of perturbations 

Let us neglect in the system of equations (84) the small contributions given by Sxy and 
Eyx- Then, for example, in the coUisionless case, we obtain for the second wave o;^ ^ klc\, 

k 

Eiyk = -j-Eizk, (93) 
Elzk — —Eixk <S Eixk- 

Thus, the second wave has a polarization of the electric field mainly along the x-axis. In 
spite of that the component Eizk ^ Ei^k, it is multiplied by a large coefficient in the first 
equation of the system (84). As a result, the contribution of this term is the same on the 
order of magnitude as that of the first term. 

In the coUisional case, the component Eizk is also defined by (93). However, its 
contribution to the first equation of the system (84) can be neglected. However, the 
contribution of the coUisional term connected with the longitudinal current in (74) which is 
proportional to Ei^k in (76) is important. 



8. DISCUSSION 

In the MHD analysis of the buoyancy instabilities, one assumes that puk/Pio <^ nuk/riiQ 
(see, e.g., Balbus 2000; Quataert 2008; Ren et al. 2009). This relation is correct for 
internal gravity waves in the neutral medium (e.g., Nekrasov 1994). It is also correct for 
perturbations rink and puk connected with the transverse perturbations Quyk 

Pilk /^ilk I r\ \ S'i i 
/ (~ Qilyk) ~ TT^ < 1 

(see 41 and 42). However, it is followed from the last equations that due to the longitudinal 
electric field perturbation E^zk which is of the order of E^zk ~ {kyCyuUdH) E^^k (see 77, 
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82, and 93) the relative pressure and number density perturbations are of the same order of 
magnitude 

Pilk jriiik I TP \ 

PiG n-io 

The ideal MHD does not capture the field Ei^^. Therefore, results obtained in the MHD 
framework and multicomponent plasma approach are different. In Sect. 4.3, we also have 
shown that V • va ~ (s'i/c^Jt'iiz- This relation is also true for internal gravity waves 
(Nekrasov 1994). The fact that V • Va 7^ for gaseous media is taken into account when 
deriving an internal energy equation (V ■ Vji is excluded from the number density and 
pressure or temperature equations) . Then in the MHD framework, one can use the condition 
of incompressibility V • Vji = in the momentum and magnetic induction equations for 
perturbations much slower than the sound waves. In our case, the divergence of the velocity 
defined by the main terms in Vjifc ~ Qnyk (see 38 and 67) and Veifc ~ Feixk (see 52 and 71) 
is also equal to zero. However, these main terms are the same for ions and electrons and 
canceled with each other at calculation of the current. Therefore, together with velocity 
perturbations proportional to the longitudinal force Fi^eizk, we must take into account 
contribution of additional small velocities connected with transverse perturbations Qiix,yk 
and Feixk- 

Prom the dispersion relation (85), we see the necessity of involving the contribution 
of values e^z, ^zx, a-nd Szz in the coUisionless case (86) (values e^z and Szx give the last 
term on the right hand-side of 88). This means that contribution of currents ji^ ~ and 
jiz ~ Eix, Eiz must be taken into account. As for the coUisional case (92), the electric field 
Eiz is not important. In the current jixki we must take into consideration the contribution 
of the current jizk as a result of collisions which is proportional to Ei^k (see 74 and 76). 
This coUisional case also is not captured by the ideal MHD. 

Thus, the standard MHD equations with simplified assumptions are not applicable 
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for the correct theory of buoyancy instabihties. Such a theory can only be given by the 
multicomponent approach used in this paper. 

The results following from (90) show that the thermal flux stabilizes the buoyancy 
instability in the case of the geometry under consideration. The instability only is possible 
in the narrow region of the temperature gradient (see Sect. 7.1). The presence of the 
background electron thermal flux (the term ~ T*q) does not play an essential role. An 
instability also is possible, if the temperature gradients of ions and electrons have the 
opposite signs. 

The contribution of collisions between electrons and ions depends on the parameter dz 
defined by (83). In the both limits (86) {dz <^ 1) and (92) {dz ^ 1), the dispersion relation 
has the same form. 

In our analysis, we have considered for generality that ions and electrons have different 
temperatures. However in (3), (6), and (7), the terms describing the energy exchange 
between species due to their collisions has not been taken into account. This is possible, if 
the dynamical timescale is smaller than the timescale of smoothing of the ion and electron 
temperatures, i.e. i/je <C fl. In the opposite case, u^^ ^, the perturbed temperatures of 
electrons and ions are almost equal one another. Equations (6) and (7) for electrons will 
keep their form because Vgi ~ Vji. In the case Teo ~ 7^0, these equations will stay the same 
with the heat flux two times less than the former one. 

Conditions of our consideration (33) and (35) are satisfied when 1 ^ Pi/H and 
1 3> kzHkypf, where pi is the ion Larmor radius. For estimations, we take u ^ g/cg 
(Teo ~ Tjo). It is obvious that these inequalities can be justified. It is easy to verify that 
conditions (79) are also satisfied. However, conditions (57) can impose some restrictions. In 
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the case Teo ~ T^o, they can be written in the form 



1 » ^ » r^, (94) 



where we have used expression for Xeo (Braginskii 1965). Prom inequahties (94), it is in 
particular followed that the case of consideration is justified when ^ ui and <^ 1. In 

the case T^q 3> T^o, the value dz is in the limits 

TeO , T„Q 1 



and can be both < 1 and > 1. 

We will further compare results obtained in this paper with the case when stratification 
is perpendicular to the magnetic field (Nekrasov and Shadmehri, 2010). However, first of 
all, we would like to say a few words about the Schwarzschild criterion of the buoyancy 
instability. It is generally accepted that this instability is possible, if the entropy increases 
in the direction of gravity. Prom a formal point of view, it is correct, if we take the 
Brunt- Vaisala frequency in the form (e.g. Balbus 2000), 

^2 ^ 1 dp d In pp-^ 

7P dz dz 

However, this expression can easily be transformed into expression (32). Thus, we see that 
the buoyancy instability exists, if the temperature increases along the gravity and the 
temperature gradient exceeds a certain threshold. 

When the thermal conduction is the dominant process in the electron temperature 
evolution, the buoyancy instability in the case g J- Bq can arise according to criterion which 
is analogous to the Schwarzschild criterion (see 65 in Nekrasov and Shadmehri, 2010). The 
same is also true when the thermal conduction is negligible (66 in Nekrasov and Shadmehri, 
2010). Both criteria are similar. Thus, we can conclude that from the point of view of 
observations it is difficult to define the role of thermal conduction in generation of instability 
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and turbulence. However, in the case g \\ Bq, the situation is different. For the negligible 
electron heat ffux, we have the similar criterion of instability as that for g ± _Bo(91). When 
the thermal conduction is dominant, the possible instability at 9 (Tjo + Tgo) /dz < is 
stabihzed (90). These results could be used by observers for determination of the mutual 
orientation of the magnetic field and gravity in astrophysical objects, e.g., in galaxy clusters. 

The true geometry of the magnetic field hues in the ICM is poorly understood. 
However, one may consider two extreme cases for the direction of the dominant magnetic 
field lines depending on the direction of gravity. The direction of the latter can have a vital 
role in driving turbulence via the possible effect of convective heat flux. Since this flux is 
mainly along the magnetic field lines, the two extreme cases are considered as either the 
magnetic field is perpendicular to the gravity or parallel to it. Thus, the true response of 
ICM to small perturbations would be possibly between these two cases. Of course, when 
the system evolves and enters into the nonlinear regime, one may expect saturation of the 
instability by rearranging the magnetic field lines. Measurements of the magnetic field in 
ICM are not consistent and there is a factor of four to ten of discrepancy depending on 
the method (e.g., Carilli and Taylor 2002). Physical mechanisms that may affect these 
observational measurements of the magnetic field in ICM are very important in this regard. 
Our analysis gives a better understanding of such a mechanism, i.e. buoyancy instability, 
though more detailed work is needed in future. 



9. CONCLUSION 



In this paper, we have investigated buoyancy instabilities in magnetized electron-ion 
plasmas with the anisotropic electron thermal flux, using the multicomponent approach 
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when the dynamical equations for the ions and electrons are solved separately via 
electric field perturbations. We have included the background electron heat fiux and 
collisions between electrons and ions. The important role of the longitudinal electric field 
perturbation, which is not captured by the ideal MHD equations, has been demonstrated. 
We have shown that the previous MHD result for the growth rate in the geometry 
considered in this paper when all background quantities are directed along the one axis is 
questionable. The reason of this has been shown to be in simplified assumptions made in 
the MHD analysis of the buoyancy instabilities and some shortcomings of the MHD. 

At the consideration of the electron heat fiux, we have adopted that the electron 
cyclotron frequency is much larger than the electron collision frequency that is typical 
for tokamaks, solar corona, and astrophysical objects such as ICM and galaxy clusters. 
The dispersion relation obtained shows that the anisotropic electron heat fiux including 
the background one stabilizes the unstable stratification except the narrow region of the 
temperature gradient. However, when gradients of the ion and electron temperatures have 
opposite signs, the medium becomes unstable. 

Results obtained in this paper are applicable to the magnetized stratified objects 
and can be useful for searching sources of turbulent transport of energy and matter. For 
astrophysical plasmas, it has been suggested that the buoyancy instability can act as a 
driving mechanism to generate turbulence in ICM and this extra source of heating may help 
to resolve the cooling flow problem (e.g., Allen 2000). However, all previous analytical and 
numerical studies are restricted to the MHD approach. Our study shows that when the true 
multifiuid nature of the system with the electron heat fiux is considered, one can not expect 
the buoyancy instability unless for a very limited range of the gradient of the temperature 
or when the gradients of the temperature of the electrons and ions have opposite signs. 
Both cases are very unlikely. However, in the case when the heat fiux does not play the 
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role, the system can be unstable due to the convective instabihty. 
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